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Abstract 

Radiative  solutions  of  the  three-dimensional  Helmholtz 

2 
equation  Au  +  k  u  =  0  are  known  to  be  representable  in  the  form 

u  =  r   exp(ikr)f(l/r,9,0)  with  f(l/r,9,0)  =  fQ(©)  +  J]     fn(9,0)r~n, 

n=l 

where  the  series  is  both  asymptotic  and  convergent.   In  fact 

f  (l/r,9,0)  is  known  to  be  an  analytic  function  of  l/r  for 

|r|  >  a  provided  all  the  singularities  of  u  are  contained  in  the 

sphere  of  radius  a.   On  the  other  hand,  in  two  dimensions,  the 

corresponding  farfield  expansion 


_,  -1/2  ikr 
u  3  r  '  e 


f  (Q)  +T]  f  (9)r_n 
n=l 


is  asymptotic,  but  is  not  convergent.  We  remedy  this  defect  by 
showing  that  any  radiative  solution  of  the  two-dimensional  Helmholtz 
equations  can  be  written  in  the  form  u  =  H^   (kr)  E(— ,9)  +  &1    ' (kr)o(-,9) , 

where  the  functions  E  and  0  are  analytic  functions  of  — ,  for  |r|  >  a, 
and  are  of  course  expressible  as  convergent  power  series  in  — .   Here 
a  is  the  radius  of  the  smallest  circle  containing  the  singularities  of 
u.   Furthermore,  we  show  how  to  deduce  E(0,Q)  and  0(0,9)  from  the 
customary  farfield  datum,  namely  f  (9) ,  and  vice  versa. 

The  convergence  of  the  new  formula  may  be  useful  in  converting 
indoor  measurements  into  farfield  patterns. 
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A  CONVERGENT  'FARFIELD'  EXPANSION  FOR  TWO-DIMENSIONAL  RADIATION  FUNCTIONS 

page  5,  line  1  -  Replace   H  (kn)   with   H  (kr) 

page  5,  equation  (7)  -  Replace   n  >  a.,       with   r  >  a, 

page  5,  line  immediately  following  equation  (8)  -  Replace   r  >  b  >  0, 
with   r  >  b  >  a  >  0, 

page  6,   equation  (ik)    -  Equation  should  read: 


(Ik)     E(l/kb  expia,   ©)   =  \ 


n=0  -nr(fexp  iaJn    Hn(kb) 


Pn(l/kbe'ia) 


page  7,  equation  (l6)    -  Replace       =  J  ,(kbc      )     \rith         =  J_-,(kbe     ) 


and     replace        =  Jo(kbcla)       with         =  JQ(kbe     ) 
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I.  INTRODUCTION 

It  is  a  known  property L  J  of  the  Helmholtz  equation, 

?     2 
V  u  +  k  u  =  0,  that  if  a  solution  u(x,y,z)  obeys  the  Sommerfeld 

radiation  condition  and  is  regular  in  the  exterior  of  a  sphere 

of  radius  a,  then  it  can  be  expanded  in  the  form 


ikr  5t^-i  f  (9,0) 
(1)       u  =  ^—  Y  A~-      :   r>a  , 


£o   r 


where  r,  9,  0  are  polar  coordinates.   The  series  is  actually 
convergent,  and  represents'  an  analytic  function  of  (l/r).   The 
functions  f  can  be  obtained  recursively1-  J  if  the  complex  radiation 
pattern  f  (9,0)  is  given. 

On  the  other  hand,  a  corresponding  theorem  in  two  dimensions 
is  unknown.   To  be  sure,  an  expansion  analogous  to  that  of  (l)  may 
be  constructed1-  J  in  the  form 


i(*r-|)y   fn(9) 
n=0 


(2)       u(x^)^V-£-  e         >     -^-      :  r  >  a  . 


r 


But  the  example  of  the  Hankel  function  H   (kr) ,  which  is  an  exact 
solution  of  the  two-dimensional  equation,  shows  that  a  series  of 
type  (2)  diverges.   An  attempt  to  repair  this  defect  by  assuming 


o.(8) 
(3)       u  ~Hw(kr)   >    -~-   :  r  >  a  , 


o       /    ,     n 

r 
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must  also  fail,  since  the  wave  function  H.J   (kr)  exp(i9)  cannot  be 
so  expressed.   Thus,  the  two-dimensional  theory  lags  behind  the 
three-dimensional  theory,  at  present.   The  purpose  of  this  paper 
is  to  remedy  this  defect. 

The  principal  result  demonstrated  below  is  the  following: 

Theorem  I;  Let  u(x,y)  =  u(r,9)  be  an  outgoing  wave  function  regular 
for  r  >  a  >  0.   Then  there  exists  a  convergent  expansion,  valid 
for  r  >  a,  in  the  form 


m    ^-,  F  (o)    (1)         ^  Gn  (e) 

(10       u  =  H^(kr)^   _£_-+H[1)(kr)  ^   -V  ' 

o    r  o    r 


where  the  series  converge  uniformly  and  absolutely  for  r  >  b  >  a, 
and  may  be  differentiated  termwise  with  respect  to  r  as  often  as 
desired. 

We  also  prove 

00  ■» 

Theorem  II:   Let  E(l/r,9)  =  \       r'^jQ)  ,  0(l/r,9)  a  y       r"nGn(9) 

0  0 

in  the  notation  of  Theorem  I.   Then  E  and  0  are  analytic  functions 
of  (l/r),  and  are  regular  in  the  interior  of  the  circle  |l/r|  <  l/b  <  l/a 
in  the  complex  (l/r)  plane,  where  b  is  any  positive  number  greater 
than  a.   We  have 

(!)  f-u*\      Wi  l-r.  a\    4.  wW. 


(k&)  u  =  HUJ(kr)   E(l/r,9)  +  H^j(kr)-0(l/r,9) 
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In  addition,  we  show  how  the  function  F  (9)  and  GQ(©)  are 
related  to  f  (9),  the  complex  radiation  pattern.   In  fact,  we  prove 

Theorem  Ills  The  leading  terms  of  E,  0,  are  related  to  the  complex 
radiation  pattern  by  the  formulae 


W      Fo(9)  =  \ 


fo(e)  +  fo(G  +  rt)l 


1Go(9)  "  \   [fo^0)  "  fo(9  +  *] 


Sinoe  the  functions  F  (9) ,  G  (0)  may  be  obtained  recursively,  as 
explained  later,  in  terms  of  F  and  G  ,  Theorem  III  makes  it  possible 
to  construct  the  convergent  expansion  (k)    in  terms  of  the  experimentally 
d  ianlngful  leading  term,  f  (9),  of  the  asymptotic  expansion  (2). 

It  is  natural  to  ask  why  the  two-dimensional  expansion  theorem 
should  turn  out  to  have  a  more  complicated  form  than  does  the 
three-dimensional  one.   The  answer  is  that  the  situation  in  two 
dimensions  is  the  more  characteristic  one.   In  fact,  it  can  be  shown 
that  in  three  dimensions  we  have  the  convergent  expansion 


(5)         u  ,  *lH  y    V9^   _  fi^\     y     Gn(9^} 

This  expansion  is  completely  analogous  to  (k)  ,  since  [_(exp  ikr)/rj 

and  H  (kr)  are  the  fundamental  solutions  in  three-  and  two -dimensions 
ox 

and  since  H.J   (kr)  is  the  derivative  of  H  (kr).   The  expansion  (5) 


can  also  be  rewritten  in  the  form  (l)  due  to  the  fact  that 
[(exp  ikr)/r]  '  =  (exp  ikr)  [ik  -  l/r)]  . 

The  relationships  discussed  in  the  last  paragraph  can  be 
made  understandable  when  the  respective  methods  of  proof  are 
explained.   Consider  first  the  proof  of  (l) .   That  proof*is  based 
on  the  expression  of  u  by  means  of  Green's  Theorem,  using  the 
fundamental  function  ((exp  ikR) /R)  and  its  normal  derivative  on 
the  sphere  r  =  a,  where  R  is  the  distance  between  the  field  point 
P(r,Q,0)  and  the  point  P(a,9  ,0  ).   Here  P  is  a  point  on  the 
sphere  r  =  a.   This  yields  an  expression  of  u  in  terms  of  integrals 
involving  u  and  its  normal  derivative  on  the  boundary.   The  elementary 
functions  of  R  which  occur  in  the  integrand  can  then  be  expanded  in 
inverse  powers  of  r,  to  prove  the  theorem.  [For  example,  we  have 


exp  ikR/R  = 


(exp  ikr) 


r  exp  ik  (R  -  r) 


and  the  second  factor  is  an  analytic  function  of  —  in  the 
neighborhood  of  infinity.]   The  proof  of  (l)  thus  depends  on  the 
elementary  character  of  the  exponential  function.   This  method  doesn't 
seem  to  be  applicable  in  two  dimensions.   The  method  of  proof  we 
have  employed  in  the  present  work  is  quite  different.   It  is  based 
on  the  fact  that  the  function  u(  x,y  )  can  be  written  in  the  form 

u  =  J>       H  (kr)  (a  cos  nb  +  bn  sin  n  9)  , 
for  r  >  a,  and,  in  addition,  on  the  well  known  fact  that  the  recurrence 
f.  [1] 


c 
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formula  for  cylinder  functions    allows  us  to  express  H  (kn)  as  a 
linear  combination  of  H  (kr)  and  H, (kr)  with,  coefficients  which  are 
polynomials  in  (l/r).   The  polynomials  are  Lommel's  polynomials  ^ -I  . 
We  show  below  that  the  resulting  expression,  namely 


(6)       u  =  H  (kr)    >   P  f-i  |  (a  cosnQ  +  b  sin  n9) 
o       /   ,   n  \kr  J         to.  n 


+  H  (kr)  \         Q    L-  )  (a  cos  n9  +  b  sin  n9)   j  r  >  b  >  a, 


can  be  rearranged  so  as  to  take  the  form  given  in  (k) . 

II.  ANALYSIS 
Suppose 

(7)       u^  +  Uyy  +  k2u  =  0  ,         n  >  a,   k  real, 


(8)       Um   \^  (tL.  -  iku)  =  0 


The  function  u  is  analytic,  for  r  >  b  >  0,  and  hence  we  may  write 

(9)       u  =   )>   — rr-T [u  (b)  cos  (n9  +  y  (b)l ,  r  >  b  >  a  >  0  , 

^  HUJ(kb)    n  n 

n=0   n  v 


(10)      u(b  cos  9,  b  sin  9)   =  u(b,9)  =  \       u  (b)  cos(n9  +  y   )    , 

n=0 


where  the  Fourier  coefficients  u  (b)  cos  y    ,   u  (b)  sin  y     tend  to  zero 

n       ' u'     nv       'n 

more  rapidly  than  any  inverse  power  of  n.   On  the  other  hand,  we  know 
that  £k~] 


(11)       Hn(kr)  =  Ho(kr)  Pn  (£)  +  H^kr)  (jjg 


where  P  (  t—  ]  ,  Q  (  t—  ]  are  the  Lommel  polynomials  |  It  J  of  degree  n 
in  the  variable  (l/kr).   In  the  notation  of  Watson,  we  have 


(12)      P  (-)    =  R   (r) 


s>(S  -  wr> 


Consider  now  the  series 


CO 


u  (b)  cos  (nQ  +  7  ) 

(i5a)        e(w,9)    =    >    -^rn^ -p>) 


H  (kb) 
n 


-<- — i  u  (b)  cos  (n9  +  7  ) 
(13b)     Q(v,G)   -   \    n  H  (Kb) ~  Q^Cv) 


0 


where  w  is  a  complex  variable.   Replace  w  by  (l/kb  exp  ia)  where 

0  <  a  <  2«;  and  b  is  a  fixed  positive  number  greater  than  a.   Then,  for 

example,  we  can  rewrite  (13&)  as  follows t 


u  (b)cos(n9  +  7)   (^(f ) 
(Ik)  E(l/kbexp  iaO)  = 


n 

ina 
e 


n=0     n  /kb 


n7V.~  exp  ia)  Hn(kb) 


-  7  - 
Now  Hurwitz  has  shown L  J  that 

(«)     m"f„  77- — -  -  '„« 

r  (v+nH-l) 


Let  z  =  kb   exp  ia  =   l/w.      Hence,  by   (12)   we  have 

la  n_1 

nP  (i/kbeia)-(^-) 
,.  ,-v               Llm             nw                      2  _     /. .    icu 

(16)  n^m      =  J„1(khc     ) 

nl 

.        n-1 

T.          n&(l/kbela)    £&§-)  ,        ic, 

Lim         XJ  2 =  j   (j^c10-) 

n  ->  00  ,  o 

n: 


r5"i 

On  the  other  hand,  for  real  b,  we  knowL^J  that 


kb^  n 

n  V 2" 

(17)      UM H  (Kb)   -  -  § 

v  ' '     n  -*  00       ,       n  it 


and  we  recall  that 


(18)     lu  (b)coe(ne  +  7  )  e"ina'|  <  o(- 
n  n  yn 


for  arbitrary  p.   These  considerations  show  that  the  series  on  the 
right  side  of  (Ik)    is  uniformly  and  absolutely  convergent.   Now  the 
series  (13&)  is  a  series  of  polynomials  in  w,  which,  as  we  have  just 
seen,  converges  uniformly  and  absolutely  for  Jw|  =  1/kb.  Hence  that 
series  converges  for  |w|  <  l/kb  and  defines  an  analytic  function. 
Thus  the  right  side  of  the  first  equation  of  (lj5)  is  convergent 


-  8  - 


for  |w|  <  —  .   The  same  is  easily  seen  to  be  true  of  the  second 
equation  of  (13).   Therefore  using  (9)  and  (ll)  we  deduce 

Lemma  I:   If  u(3?,Q)  is  an  outgoing  wave  function  in  two  variables, 
regular  outside  a  circle  r  =  a,  then 


S— ,  u  (b)cos(a©  +  7  ) 

(19)      u  =  H  (kr)  >  -5 R   (kr) 

0     Zj     H  (a)         n'° 
n=0      n 


^  u  (b)cos(n0+7  ) 

+  H_(kr)  >  -S _2  R     (kr) 

1    Zj    H  (kb)       n~1}1 

n=0      n 


where  R  n  and  R      are  polynomials  of  degree  n  in  the  variable  (ifisr)  ,   and 
1  ■  '   n,U  ~ ~~~  n—XjJ_ 

where  the  series  converge  uniformly  and  absolutely  for  r  >  b  >  a. 
The  series  also  converge  uniformly  and  absolutely  for  a  complex  r, 
if  [r|  >  b  >  a. 

This  Lemma  easily  implies  Theorem  II  of  the  Introduction,  since 
the  uniform  convergence  of  the  series  for  E  and  0,  of  (13&)  a-ncl  (ljb) 
shows  that  E  and  0  are  analytic  functions  of  (l/r)  in  the  circle 
(l/r)  <  l/b  in  the  complex  l/r  plane.   Theorem  I  is  a  trivial 
corollary  of  Theorem  II. 

III.  DETERMINATION  OF  THE  COEFFICIENTS 
We  have  thus  far  demonstrated  the  existence  of  a  convergent 
expansion  of  the  form  given  in  equation  (k) .      It  is  now  appropriate 
to  show  how  to  determine  the  coefficients  in  this  expansion.   We 
consider  two  cases:  In  the  first  case  u  is  given  on  some  circle  r  =  b  >  a. 


In  the  second  case  the  complex  far  field  amplitude  f  (9)  is  given. 

The  first  case  is  trivial.   In  fact,  the  series  (9)  will 
supply  the  coefficient  u  (b)  and  the  phase  angles  7  .   These  can  be 
inserted  in  (19)  and  the  resulting  series  can  be  rearranged  in 
descending  powers  of  r,  by  the  Weierstrass  double-series  theorem. 

It  Is  of  greater  interest  to  construct  the  series,  (k)   when 
we  are  given  data  at  r  =  00.   Suppose  therefore  that  we  are  given  the 
complex  far  field  amplitude  f  (9),  where,  by  (2) 


(20)      J?m       ^  u  exp(-ikr  +  f )   -   f  (9) 


TJ 


Notice  that  F  (0)  =  E(0,Q)  and  G  (9)  =  0(0,9).   It  is  easily  seen 

that,  if  the  functions  F  (9),  G  (9)  of  equation  (l^)  are  given, 

then  F  and  G  can  be  computed  recursively.   To  see  this  fact, 

simply  substitute  (K)    into  the  equation  v  u+k  u  =  0,  and  express 

all  derivatives  of  H  and  EL  in  terms  of  H  and  EL  ,  by  means  of  the 

o      1  o      1 

recurrence  formulae , 


(21) 


1  11     2. 

H   =  -H,  ,      H  =  +  -  H  -  H 
o      1        o     r  1   o 


t        II 


2H     -2H       „ 

— i.  =  °_    Hl   =  -(-  H.-H 

r       r   '        1     vr  1   o' 


Then  collect  the  terms  multiplying  H  and  H  .   The  net  coefficients 

of  H  and  H  must  vanish  identically,  and  this  gives  expressions  for 
o      1 

F  ,G  in  terms  of  F  and  G  .   For  example  we  find 
n'  n  00 


-  10  - 


(22) 


Fl  =  "  |<°o"  +  Go} 


1   " 
1     2  o 


Therefore  the  problem  is  really  to  find  F  and  G  ,  given  f  (©).   In 

order  to  find  F  and  G  ,  i.e.-,  to  deduce  the  result  given  in  Theorem 
o      o'     ' 

III,  we  notice  first  that,  as  r  -+  e», 


(21a)     H(l)  ~   (-I)n  5  (1) 
v  '  '      n      v  '        o 


(21b)     H^2)  -   (+i)U  Ho(2)  , 


*(«.   ZZ  .1(kr-?}   ~,  (2):.AH,-i(kr-^ 


where 


(21c)     H  v  '   a  A/-T—  e        j  H  v  ;  =  V-^:  e 
v    '      o      V  itkr  o      V  jtkt 


Also,  we  have 

(22)  H^m)(kr)Pn  +  H^kr)^  =  H^m)  (kr)  ,   m  =  1,2. 

These  equations  show  that  as  r  ->  °° 

(23)  Pn  -  iQ^  -  (-i)n  ,  Pn  +  iQ^  ->  (+i)n  , 
and,  consequently,  that 


11 


> 


(21+) 


P„     -  (~1)P      :     Po     ■,   -  0 
2p        v      '  2p+l 


for  r  -*  oo 


%-° 


Vi  -  (-l): 


Therefore,  by  (19)  and  (21a)  and  (22),  we  have 


(25) 


fo(G) 


un(b)(-i)n 

v-,  \ cos(n©  +  y 

(1)/..  n  a' 


H^(kb) 

n=0    n 


But,  since  F  (9)  =  E(0,9),  G  (9)  =  0(0,9) ,  we  find,  by  letting 
r  -»  co  in  (ka.)    and  by  using  (21a)  and  (k)  ,   that 


(26)      fQ(9)   =   Fq(9)  -  iGQ(9)   =  E(0,9)  -  10(0,9) 


Furthermore,  if  we  use  (19)  and  (12),  and  exploit  the  simplification 
introduced  by  (2*0  ,  we  find  that  F  (9)  contains  only  even -numbered 
Fourier  coefficients,  and  G  (9)  contains  only  odd  ones.   In  fact 


(2  7)      Fq(9)   = 


upT1  (b)(-l)P 

H*  (kb)     COS  (2p9  +  W 
2p 


G  (9)   = 


V   _2P±i 


u£>   (-DP 


J    H2p+l(kb) 


cos 


[(2p+l)9  +  7o^J 


2p+l- 


If  we  compare  (25)  with  (27),  we  find,  in  view  of  (26),  the 
result  expressed  in  Theorem  III  of  the  introduction.   Thus  we  have 

expressed  F  ,G  in  terms  of  f  ,  as  we  desired. 
00  o 
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